Sampling Phase Space Dividing Surfaces Constructed from Normally Hyperbolic Invariant Manifolds (NHIMs).

In this paper, we further investigate the construction of a phase space dividing surface (DS) from a normally hyperbolic invariant manifold (NHIM) and the sampling procedure for the resulting dividing surface described in earlier work ( Wiggins , S. ; J. Chem. Phys. 2016 , 144 , 054107 ). Our discussion centers on the relationship between geometrical structures and dynamics for 2 and 3 degree of freedom (DoF) systems, specifically, the construction of a DS from a NHIM. We show that if the equation for the NHIM and associated DS is known (e.g., as obtained from Poincaré-Birkhoff normal form theory), then the numerical procedure described in Wiggins et al. ( J. Chem. Phys. 2016 , 144 , 054107 ) gives the same result as a sampling method based upon the explicit form of the NHIM. After describing the sampling procedure in a general context, it is applied to a quadratic Hamiltonian normal form near an index-one saddle where explicit formulas exist for both the NHIM and the DS. It is shown for both 2 and 3 DoF systems that a version of the general sampling procedure provides points on the analytically defined DS with the correct microcanonical density on the constant-energy DS. Excellent agreement is obtained between analytical and numerical averages of quadratic energy terms over the DS for a range of energies.

Nonadiabatic semiclassical dynamics in the mixed quantum-classical initial value representation.

We extend the Mixed Quantum-Classical Initial Value Representation (MQC-IVR), a semiclassical method for computing real-time correlation functions, to electronically nonadiabatic systems using the Meyer-Miller-Stock-Thoss (MMST) Hamiltonian in order to treat electronic and nuclear degrees of freedom (dofs) within a consistent dynamic framework. We introduce an efficient symplectic integration scheme, the MInt algorithm, for numerical time evolution of the phase space variables and monodromy matrix under the non-separable MMST Hamiltonian. We then calculate the probability of transmission through a curve crossing in model two-level systems and show that MQC-IVR reproduces quantum-limit semiclassical results in good agreement with exact quantum methods in one limit, and in the other limit yields results that are in keeping with classical limit semiclassical methods like linearized IVR. Finally, exploiting the ability of the MQC-IVR to quantize different dofs to different extents, we present a detailed study of the extents to which quantizing the nuclear and electronic dofs improves numerical convergence properties without significant loss of accuracy.

Empirical Classification of Trajectory Data: An Opportunity for the Use of Machine Learning in Molecular Dynamics.

Classical Hamiltonian trajectories are initiated at random points in phase space on a fixed energy shell of a model two degrees of freedom potential, consisting of two interacting minima in an otherwise flat energy plane of infinite extent. Below the energy of the plane, the dynamics are demonstrably chaotic. However, most of the work in this paper involves trajectories at a fixed energy that is 1% above that of the plane, in which regime the dynamics exhibit behavior characteristic of chaotic scattering. The trajectories are analyzed without reference to the potential, as if they had been generated in a typical direct molecular dynamics simulation. The questions addressed are whether one can recover useful information about the structures controlling the dynamics in phase space from the trajectory data alone, and whether, despite the at least partially chaotic nature of the dynamics, one can make statistically meaningful predictions of trajectory outcomes from initial conditions. It is found that key unstable periodic orbits, which can be identified on the analytical potential, appear by simple classification of the trajectories, and that the specific roles of these periodic orbits in controlling the dynamics are also readily discerned from the trajectory data alone. Two different approaches to predicting trajectory outcomes from initial conditions are evaluated, and it is shown that the more successful of them has ∼90% success. The results are compared with those from a simple neural network, which has higher predictive success (97%) but requires the information obtained from the "by-hand" analysis to achieve that level. Finally, the dynamics, which occur partly on the very flat region of the potential, show characteristics of the much-studied phenomenon called "roaming." On this potential, it is found that roaming trajectories are effectively "failed" periodic orbits and that angular momentum can be identified as a key controlling factor, despite the fact that it is not a strictly conserved quantity. It is also noteworthy that roaming on this potential occurs in the absence of a "roaming saddle," which has previously been hypothesized to be a necessary feature for roaming to occur.

Roaming: A Phase Space Perspective.

In this review we discuss the recently described roaming mechanism for chemical reactions from the point of view of nonlinear dynamical systems in phase space. The recognition of the roaming phenomenon shows the need for further developments in our fundamental understanding of basic reaction dynamics, as is made clear by considering some questions that cut across most studies of roaming: Is the dynamics statistical? Can transition state theory be applied to estimate roaming reaction rates? What role do saddle points on the potential energy surface play in explaining the behavior of roaming trajectories? How do we construct a dividing surface that is appropriate for describing the transformation from reactants to products for roaming trajectories? How should we define the roaming region? We show that the phase space perspective on reaction dynamics provides the setting in which these questions can be properly framed and answered. We illustrate these ideas by considering photodissociation of formaldehyde. The phase-space formulation allows an unambiguous description of all possible reactive events, which also allows us to uncover the phase space mechanism that explains which trajectories roam, as opposed to evolving toward a different reactive event.

Phase space barriers and dividing surfaces in the absence of critical points of the potential energy: Application to roaming in ozone.

We examine the phase space structures that govern reaction dynamics in the absence of critical points on the potential energy surface. We show that in the vicinity of hyperbolic invariant tori, it is possible to define phase space dividing surfaces that are analogous to the dividing surfaces governing transition from reactants to products near a critical point of the potential energy surface. We investigate the problem of capture of an atom by a diatomic molecule and show that a normally hyperbolic invariant manifold exists at large atom-diatom distances, away from any critical points on the potential. This normally hyperbolic invariant manifold is the anchor for the construction of a dividing surface in phase space, which defines the outer or loose transition state governing capture dynamics. We present an algorithm for sampling an approximate capture dividing surface, and apply our methods to the recombination of the ozone molecule. We treat both 2 and 3 degrees of freedom models with zero total angular momentum. We have located the normally hyperbolic invariant manifold from which the orbiting (outer) transition state is constructed. This forms the basis for our analysis of trajectories for ozone in general, but with particular emphasis on the roaming trajectories.

Toward Understanding the Roaming Mechanism in H + MgH → Mg + HH Reaction.

The roaming mechanism in the reaction H + MgH →Mg + HH is investigated by classical and quantum dynamics employing an accurate ab initio three-dimensional ground electronic state potential energy surface. The reaction dynamics are explored by running trajectories initialized on a four-dimensional dividing surface anchored on three-dimensional normally hyperbolic invariant manifold associated with a family of unstable orbiting periodic orbits in the entrance channel of the reaction (H + MgH). By locating periodic orbits localized in the HMgH well or involving H orbiting around the MgH diatom, and following their continuation with the total energy, regions in phase space where reactive or nonreactive trajectories may be trapped are found. In this way roaming reaction pathways are deduced in phase space. Patterns similar to periodic orbits projected into configuration space are found for the quantum bound and resonance eigenstates. Roaming is attributed to the capture of the trajectories in the neighborhood of certain periodic orbits. The complex forming trajectories in the HMgH well can either return to the radical channel or "roam" to the MgHH minimum from where the molecule may react.

Phase Space Structures Explain Hydrogen Atom Roaming in Formaldehyde Decomposition.

We re-examine the prototypical roaming reaction--hydrogen atom roaming in formaldehyde decomposition--from a phase space perspective. Specifically, we address the question "why do trajectories roam, rather than dissociate through the radical channel?" We describe and compute the phase space structures that define and control all possible reactive events for this reaction, as well as provide a dynamically exact description of the roaming region in phase space. Using these phase space constructs, we show that in the roaming region, there is an unstable periodic orbit whose stable and unstable manifolds define a conduit that both encompasses all roaming trajectories exiting the formaldehyde well and shepherds them toward the H2···CO well.

Reaction Path Bifurcation in an Electrocyclic Reaction: Ring-Opening of the Cyclopropyl Radical.

Following previous work [J. Chem. Phys. 2013, 139, 154108] on a simple model of a reaction with a post-transition state valley ridge inflection point, we study the chemically important example of the electrocyclic cyclopropyl radical ring-opening reaction using direct dynamics and a reduced dimensional potential energy surface. The overall reaction requires con- or disrotation of the methylenes, but the initial stage of the ring-opening involves substantial internal rotation of only one methylene. The reaction path bifurcation is then associated with the relative sense of rotation of the second methylene. Clear deviations of reactive trajectories from the disrotatory intrinsic reaction coordinate (IRC) for the ring-opening are observed and the dynamical mechanism is discussed. Several features observed in the model system are found to be preserved in the more complex and higher dimensional ring-opening reaction. Most notable is the sensitivity of the reaction mechanism to the shape of the potential manifested as a Newtonian kinetic isotope effect upon deuterium substitution of one of the methylene hydrogens. Dependence of the product yield on frictional dissipation representing external environmental effects is also presented. The dynamics of the post-transition state cyclopropyl radical ring-opening are discussed in detail, and the use of low dimensional models as tools to analyze complicated organic reaction mechanisms is assessed in the context of this reaction.

Nonstatistical dynamics on the caldera.

We explore both classical and quantum dynamics of a model potential exhibiting a caldera: that is, a shallow potential well with two pairs of symmetry related index one saddles associated with entrance/exit channels. Classical trajectory simulations at several different energies confirm the existence of the "dynamical matching" phenomenon originally proposed by Carpenter, where the momentum direction associated with an incoming trajectory initiated at a high energy saddle point determines to a considerable extent the outcome of the reaction (passage through the diametrically opposing exit channel). By studying a "stretched" version of the caldera model, we have uncovered a generalized dynamical matching: bundles of trajectories can reflect off a hard potential wall so as to end up exiting predominantly through the transition state opposite the reflection point. We also investigate the effects of dissipation on the classical dynamics. In addition to classical trajectory studies, we examine the dynamics of quantum wave packets on the caldera potential (stretched and unstretched). These computations reveal a quantum mechanical analogue of the "dynamical matching" phenomenon, where the initial expectation value of the momentum direction for the wave packet determines the exit channel through which most of the probability density passes to product.

Roaming dynamics in ion-molecule reactions: phase space reaction pathways and geometrical interpretation.

A model Hamiltonian for the reaction CH4(+) -> CH3(+) + H, parametrized to exhibit either early or late inner transition states, is employed to investigate the dynamical characteristics of the roaming mechanism. Tight/loose transition states and conventional/roaming reaction pathways are identified in terms of time-invariant objects in phase space. These are dividing surfaces associated with normally hyperbolic invariant manifolds (NHIMs). For systems with two degrees of freedom NHIMS are unstable periodic orbits which, in conjunction with their stable and unstable manifolds, unambiguously define the (locally) non-recrossing dividing surfaces assumed in statistical theories of reaction rates. By constructing periodic orbit continuation/bifurcation diagrams for two values of the potential function parameter corresponding to late and early transition states, respectively, and using the total energy as another parameter, we dynamically assign different regions of phase space to reactants and products as well as to conventional and roaming reaction pathways. The classical dynamics of the system are investigated by uniformly sampling trajectory initial conditions on the dividing surfaces. Trajectories are classified into four different categories: direct reactive and non-reactive trajectories, which lead to the formation of molecular and radical products respectively, and roaming reactive and non-reactive orbiting trajectories, which represent alternative pathways to form molecular and radical products. By analysing gap time distributions at several energies, we demonstrate that the phase space structure of the roaming region, which is strongly influenced by nonlinear resonances between the two degrees of freedom, results in nonexponential (nonstatistical) decay.

Nonstatistical dynamics on potentials exhibiting reaction path bifurcations and valley-ridge inflection points.

We study reaction dynamics on a model potential energy surface exhibiting post-transition state bifurcation in the vicinity of a valley ridge inflection (VRI) point. We compute fractional yields of products reached after the VRI region is traversed, both with and without dissipation. It is found that apparently minor variations in the potential lead to significant changes in the reaction dynamics. Moreover, when dissipative effects are incorporated, the product ratio depends in a complicated and highly non-monotonic fashion on the dissipation parameter. Dynamics in the vicinity of the VRI point itself play essentially no role in determining the product ratio, except in the highly dissipative regime.

Isomerization dynamics of a buckled nanobeam.

We analyze the dynamics of a model of a nanobeam under compression. The model is a two-mode truncation of the Euler-Bernoulli beam equation subject to compressive stress applied at both ends. We consider parameter regimes where the first mode is unstable and the second mode can be either stable or unstable, and the remaining modes (neglected) are always stable. Material parameters used correspond to a silicon nanobeam. The two-mode model Hamiltonian is the sum of a (diagonal) kinetic energy term and a potential energy term. The form of the potential energy function suggests an analogy with isomerization reactions in chemistry, where "isomerization" here corresponds to a transition between two stable beam configurations. We therefore study the dynamics of the buckled beam using the conceptual framework established for the theory of isomerization reactions. When the second mode is stable the potential energy surface has an index one saddle, and when the second mode is unstable the potential energy surface has an index two saddle and two index one saddles. Symmetry of the system allows us to readily construct a phase space dividing surface between the two "isomers" (buckled states); we rigorously prove that, in a specific energy range, it is a normally hyperbolic invariant manifold. The energy range is sufficiently wide that we can treat the effects of the index one and index two saddles on the isomerization dynamics in a unified fashion. We have computed reactive fluxes, mean gap times, and reactant phase space volumes for three stress values at several different energies. In all cases the phase space volume swept out by isomerizing trajectories is considerably less than the reactant density of states, proving that the dynamics is highly nonergodic. The associated gap time distributions consist of one or more "pulses" of trajectories. Computation of the reactive flux correlation function shows no sign of a plateau region; rather, the flux exhibits oscillatory decay, indicating that, for the two-mode model in the physical regime considered, a rate constant for isomerization does not exist.

Index k saddles and dividing surfaces in phase space with applications to isomerization dynamics.

In this paper, we continue our studies of the phase space geometry and dynamics associated with index k saddles (k > 1) of the potential energy surface. Using Poincaré-Birkhoff normal form (NF) theory, we give an explicit formula for a "dividing surface" in phase space, i.e., a codimension one surface (within the energy shell) through which all trajectories that "cross" the region of the index k saddle must pass. With a generic non-resonance assumption, the normal form provides k (approximate) integrals that describe the saddle dynamics in a neighborhood of the index k saddle. These integrals provide a symbolic description of all trajectories that pass through a neighborhood of the saddle. We give a parametrization of the dividing surface which is used as the basis for a numerical method to sample the dividing surface. Our techniques are applied to isomerization dynamics on a potential energy surface having four minima; two symmetry related pairs of minima are connected by low energy index 1 saddles, with the pairs themselves connected via higher energy index 1 saddles and an index 2 saddle at the origin. We compute and sample the dividing surface and show that our approach enables us to distinguish between concerted crossing ("hilltop crossing") isomerizing trajectories and those trajectories that are not concerted crossing (potentially sequentially isomerizing trajectories). We then consider the effect of additional "bath modes" on the dynamics, by a study of a four degree-of-freedom system. For this system we show that the normal form and dividing surface can be realized and sampled and that, using the approximate integrals of motion and our symbolic description of trajectories, we are able to choose initial conditions corresponding to concerted crossing isomerizing trajectories and (potentially) sequentially isomerizing trajectories.

Phase space structure and dynamics for the Hamiltonian isokinetic thermostat.

We investigate the phase space structure and dynamics of a Hamiltonian isokinetic thermostat, for which ergodic thermostat trajectories at fixed (zero) energy generate a canonical distribution in configuration space. Model potentials studied consist of a single bistable mode plus transverse harmonic modes. Interpreting the bistable mode as a reaction (isomerization) coordinate, we establish connections with the theory of unimolecular reaction rates, in particular the formulation of isomerization rates in terms of gap times. In the context of molecular reaction rates, the distribution of gap times (or associated lifetimes) for a microcanonical ensemble initiated on the dividing surface is of great dynamical significance; an exponential lifetime distribution is usually taken to be an indicator of "statistical" behavior. Moreover, comparison of the magnitude of the phase space volume swept out by reactive trajectories as they pass through the reactant region with the total phase space volume (classical density of states) for the reactant region provides a necessary condition for ergodic dynamics. We compute gap times, associated lifetime distributions, mean gap times, reactive fluxes, reactive volumes, and total reactant phase space volumes for model thermostat systems with three and four degrees of freedom at three different temperatures. At all three temperatures, the necessary condition for ergodicity is approximately satisfied. At high temperatures a nonexponential lifetime distribution is found, while at low temperatures the lifetime is more nearly exponential. The degree of exponentiality of the lifetime distribution is quantified by computing the information entropy deficit with respect to pure exponential decay. The efficacy of the Hamiltonian isokinetic thermostat is examined by computing coordinate distributions averaged over single long trajectories initiated on the dividing surface.

Bulgac-Kusnezov-Nosé-Hoover thermostats.

In this paper, we formulate Bulgac-Kusnezov constant temperature dynamics in phase space by means of non-Hamiltonian brackets. Two generalized versions of the dynamics are similarly defined, one where the Bulgac-Kusnezov demons are globally controlled by means of a single additional Nosé variable, and another where each demon is coupled to an independent Nosé-Hoover thermostat. Numerically stable and efficient measure-preserving time-reversible algorithms are derived in a systematic way for each case. The chaotic properties of the different phase space flows are numerically illustrated through the paradigmatic example of the one-dimensional harmonic oscillator. It is found that, while the simple Bulgac-Kusnezov thermostat is apparently not ergodic, both of the Nosé-Hoover controlled dynamics sample the canonical distribution correctly.

Microcanonical rates, gap times, and phase space dividing surfaces.

The general approach to classical unimolecular reaction rates due to Thiele is revisited in light of recent advances in the phase space formulation of transition state theory for multidimensional systems. Key concepts, such as the phase space dividing surface separating reactants from products, the average gap time, and the volume of phase space associated with reactive trajectories, are both rigorously defined and readily computed within the phase space approach. We analyze in detail the gap time distribution and associated reactant lifetime distribution for the isomerization reaction HCN <==> CNH, previously studied using the methods of phase space transition state theory. Both algebraic (power law) and exponential decay regimes have been identified. Statistical estimates of the isomerization rate are compared with the numerically determined decay rate. Correcting the RRKM estimate to account for the measure of the reactant phase space region occupied by trapped trajectories results in a drastic overestimate of the isomerization rate. Compensating but as yet not fully understood trapping mechanisms in the reactant region serve to slow the escape rate sufficiently that the uncorrected RRKM estimate turns out to be reasonably accurate, at least at the particular energy studied. Examination of the decay properties of subensembles of trajectories that exit the HCN well through either of two available symmetry related product channels shows that the complete trajectory ensemble effectively attains the full symmetry of the system phase space on a short time scale t approximately < 0.5 ps, after which the product branching ratio is 1:1, the "statistical" value. At intermediate times, this statistical product ratio is accompanied by nonexponential (nonstatistical) decay. We point out close parallels between the dynamical behavior inferred from the gap time distribution for HCN and nonstatistical behavior recently identified in reactions of some organic molecules.

Reversible measure-preserving integrators for non-Hamiltonian systems.

We present a systematic method for deriving reversible measure-preserving integrators for non-Hamiltonian systems such as the Nosé-Hoover thermostat and generalized Gaussian moment thermostat (GGMT). Our approach exploits the (non-Poisson) bracket structure underlying the thermostat equations of motion. Numerical implementation for the GGMT system shows that our algorithm accurately conserves the thermostat energy function. We also study position and momentum distribution functions obtained using our integrator.

Classical and quantum mechanics of diatomic molecules in tilted fields.

We investigate the classical and quantum mechanics of diatomic molecules in noncollinear (tilted) static electric and nonresonant linearly polarized laser fields. The classical diatomic in tilted fields is a nonintegrable system, and we study the phase space structure for physically relevant parameter regimes for the molecule KCl. While exhibiting low-energy (pendular) and high-energy (free-rotor) integrable limits, the rotor in tilted fields shows chaotic dynamics at intermediate energies, and the degree of classical chaos can be tuned by changing the tilt angle. We examine the quantum mechanics of rotors in tilted fields. Energy-level correlation diagrams are computed, and the presence of avoided crossings quantified by the study of nearest-neighbor spacing distributions as a function of energy and tilting angle. Finally, we examine the influence of classical periodic orbits on rotor wave functions. Many wave functions in the tilted field case are found to be highly nonseparable in spherical polar coordinates. Localization of wave functions in the vicinity of classical periodic orbits, both stable and unstable, is observed for many states.

Semiclassical calculation of the vibrational echo.

The infrared echo measurement probes the time scales of the molecular motions that couple to a vibrational transition. Computation of the echo observable within rigorous quantum mechanics is problematic for systems with many degrees of freedom, motivating the development of semiclassical approximations to the nonlinear optical response. We present a semiclassical approximation to the echo observable, based on the Herman-Kluk propagator. This calculation requires averaging over a quantity generated by two pairs of classical trajectories and associated stability matrices, connected by a pair of phase-space jumps. Quantum, classical, and semiclassical echo calculations are compared for a thermal ensemble of noninteracting anharmonic oscillators. The semiclassical approach uses input from classical mechanics to reproduce the significant features of a complete, quantum mechanical calculation of the nonlinear response.

Global analysis of periodic orbit bifurcations in coupled Morse oscillator systems: time-reversal symmetry, permutational representations and codimension-2 collisions.

In this paper we study periodic orbit bifurcation sequences in a system of two coupled Morse oscillators. Time-reversal symmetry is exploited to determine periodic orbits by iteration of symmetry lines. The permutational representation of Tsuchiya and Jaffe is employed to analyze periodic orbit configurations on the symmetry lines. Local pruning rules are formulated, and a global analysis of possible bifurcation sequences of symmetric periodic orbits is made. Analysis of periodic orbit bifurcations on symmetry lines determines bifurcation sequences, together with periodic orbit periodicities and stabilities. The correlation between certain bifurcations is explained. The passage from an integrable limit to nointegrability is marked by the appearance of tangent bifurcations; our global analysis reveals the origin of these ubiquitous tangencies. For period-1 orbits, tangencies appear by a simple disconnection mechanism. For higher period orbits, a different mechanism involving 2-parameter collisions of bifurcations is found. (c) 1999 American Institute of Physics.